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Cosmic ray transport in MHD turbulence 

Huirong Yan ^ and A. Lazarian ^ 

Abstract. Numerical simulations shed light onto earlier not trackable prob- 
lem of magnetohydrodynamic (MHD) turbulence. They allowed to test the pre- 
dictions of different models and choose the correct ones. Inevitably, this progress 
calls for revisions in the picture of cosmic ray (CR) transport. It also shed light 
on the problems with the present day numerical modeling of CR. In this paper 
we focus on the analytical way of describing CR propagation and scattering, 
which should be used in synergy with the numerical studies. In particular, we 
use recently established scaling laws for MHD modes to obtain the transport 
properties for CRs. We include nonlinear effects arising from large scale trap- 
ping, to remove the 90 degree divergence. We determine how the efficiency of the 
scattering and CR mean free path depend on the characteristics of ionized me- 
dia, e.g. plasma /?, Coulomb coUisional mean free path. Implications for particle 
transport in interstellar medium and solar corona are discussed. We also exam- 
ine the perpendicular transport of CRs. Perpendicular transport depends on the 
comparison of parallel mean free path and the injection scale of the turbulence, 
as well as the Alfvenic Mach number. Normal turbulence does not allow subd- 
iffusion unless there are slab waves. The critical scale below which subdiffusion 
applies is provided. These results can be used to compare with the numerical 
calculations, provided that these calculations use the structure of magnetic field 
which is consistent with the numerical studies of MHD turbulence. 



1. Introduction 

The propagation and acceleration of cosmic rays (CRs) is governed by their 
interactions with magnetic fields. The perturbations of turbulent magnetic field 
can be accounted for by direct numerical scattering simulations (see Giacalone 
&: Jokipii 1999) or by quasi-linear theory, QLT (see Jokipii 1966, Schlickeiser 
2002). The problem with direct numerical simulations of scattering is that the 
present-day MHD simulations have rather limited inertial range. As a result, 
they are sufficient to test theoretical expectations (see Cho &: Lazarian 2003) of 
the power spectra and its anisotropy, but are not adequate to study the particle 
scattering. Indeed, if one directly uses the results of MHD simulations, it is 
likely that the scattering will be happening by magnetic perturbations beyond 
the inertial range. 

Alternative way would be creating magnetic fields in synthetic data cubes. 
However, both theory (Goldreich &: Sridhar 1995, henceforth GS95) and numer- 
ical simulations (Cho k, Vishniac 2000, Maron & Goldreich 2001, Cho, Lazarian 
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k Vishniac 2002, Muller & Biskamp 2000, Cho k Lazarian 2003 and ref. therein) 
show that the Alfvenic modes exhibit scale-dependent anisotropy, i.e. the de- 
gree of anisotropy of eddies changes with scale. The complication here is that 
this anisotropy can be seen only in the system of reference connected with the 
local magnetic field, which makes the accepted procedures of generating random 
fields in Fourier space inapplicable. As CRs do sample local magnetic fields, 
it seems essential to use the correct description of this field, which presents an 
exciting, but still unsolved problem for numerical modeling of CR propagation. 
In this situation, analytical calculations are indispensable. These studies should 
also provide one with guidance in interpreting the present-day CR propagation 
simulations and designing the new ones. 

So far the most frequently used analytical formalism is QLT. While QLT 
provides simple physical insights into scattering, it is known to have problems. 
For instance, it fails in treating 90° scattering and perpendicular transport. 
Indeed, many attempts have been made to improve the QLT and various non- 
linear theories have been attempted (see Dupree 1966, Volk 1973, Jones, Kaiser 
& Birmingham 1973, Goldstein 1976, Matthaeus et al. 2003, Shalchi 2005). Most 
of the analysis so far are confined to slab model and 2D of MHD perturbations. 
We feel that it is important to extend the work on non-linear treatment of CR 
scattering to models of large scale MHD turbulence. Here we use models that arc 
motivated by the recent studies of MHD turbulence (GS95, see Cho, Lazarian, 
Vishniac 2003 for a review and references therein). We showed in Yan & lazarian 
(2002, 2004, henceforth YL02, 04) that scattering of CRs is dominated by fast 
modes. We shall reexamine our earlier conclusions. More important, by solving 
the 90° scattering, we shall determine how CR mean free path varies in different 
environments. 

We shall also address the problem of perpendicular transport of CR, which is 
another important problem in which QLT encounters serious difficulties. Indeed, 
the idea of CR transport in the direction perpendicular to the mean magnetic 
field being dominated by the field line random walk (FLRW, Jokipii 1966, Jokipii 
& Parker 1969, Forman et al. 1974) only applies to a restricted situation where 
the turbulence perturbations are small and CRs' mean free paths are larger than 
the injection scale of MHD turbulence. It was speculated that in some regime the 
particle motions can be subdiffusive (Kota &l Jokipii 2000, Mace et al 2000, Qin 
at al. 2002a, Shalchi 2005, Webb et al. 2006). This could indicate a substantial 
shift in the paradigm of CR transport. How realistic is the subdiffusion in the 
presence of turbulence injected on large scales? In this paper we shall use the 
tested model of turbulence to study the issue of subdiffusion. 



2. Nonlinear transport in MHD turbulence 

The basic assumption of quasi-linear theory is that particles follow unperturbed 
orbits until they are scattered off by ^ 90°. In reality, particle's pitch an- 
gle varies gradually with the variation of magnetic field due to conservation of 
adiabatic invariant. The average uncertainty of parallel speed Av[| is Av[[ ~ 

^< > / Bq^^ , where Bq is mean magnetic field strength. Perturbation 
5B\\ exists owing to slow modes in incompressible turbulence and high B = □''"^ 

II ^maq 
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turbulence. In low /3 compressible turbulence, 5B^^ arises from fast modes. This 
results in broadening of the resonance. Because of the dispersion of the pitch 
angle Afi and therefore parallel speed Auy, the guiding center is perturbed about 
the mean position < z >= vfit as they move along the field lines. 




Figure 1. Scattering of CRs in strong turbulence (Alfvenic Mach number 
Ma = 6B^/Bq ~ 1). Left: scattering coefficient in strong incompressible 
turbulence; right: scattering coefficient in fast mode turbulence. 

The characteristic phase function e*^ll^^*^ deviates from that for plane waves, 
assuming the guiding center has a Gaussian distribution along the field line, we 
get the broadened resonance function 



2tt 



27r|A;i|At;ii 



■ exp 



2(A;|| v/i — u + nQ)^ 



(1) 



in contrast to the S{k\\v\\ — uj zizQ) function resonance in QLT. In the equation, 
A;|| is the component of wave vector k projected along the local magnetic field, 
CO is the wave frequency and is the Larmor frequency of the particle. The 
corresponding pitch angle diffusion can then be obtained by replacing the 6 
function with the new resonance function. For gyroresonance (n = ±1,2,...), 
the result is similar to that from QLT. This is reasonable as the gyroresonance 
works on a local scale so that the effect of large scale trapping is negligible. 

On the other hand, the dispersion of the t^y means CRs with a much wider 
range of pitch angle can be scattered through TTD (transit time damping, n=0) 
by the compressible modes. In incompressible medium, TTD happens with 
pseudo-Alfven modes (the incompressible limit of slow modes). The scattering 
by TTD produces a flat curve at large pitch angle (including 90°, see Fig(T]), 
which is conflrmed by our preliminary result from simulations (Beresnyak et al. 
in preparation). 

Recent studies have shown that Alfven (and slow) modes exhibit scale- 
dependent anisotropy and follow GS95 relations. The mixing motions associated 

with Alfvenic turbulence induce the scale-dependent anisotropy /cy ~ k'^^L'^^^ 
on slow modes, which on their own would evolve on substantially longer time 
scale. Fast modes in low Pp plasma, on the other hand, develop on their own. 
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Figure 2. Left: The turbulence truncation scales in Galactic halo and warm 
ionized medium (WIM). The damping curve becomes flat around 90° in halo 
due to field line wandering; For WIM, both viscous and coUisionless damping 
are applicable, right: The mean free paths in halo (solid line) and WIM 
(dashed line). 



as their phase velocity is only marginally affected by mixing motions induced 
by Alfven modes. According to Cho & Lazarian (2002), fast modes follow an 

3 

isotropic "acoustic" cascade with one dimensional energy spectrum W{k) oc k~2 . 
We adopt these models of MHD turbulence in the paper. 

Although the 90° divergence can be removed in the nonlinear treatment, 
the mean free path would be still unrealistically large if only Alfven and slow 
modes are counted. This is because at small pitch angle (propagating nearly 
parallel to the magnetic field), only gyroresonance applies, which is dominated 
by the interaction with fast modes (YL02, 04). 

In reality, compressible turbulence is subjected to damping. The damping 
depends on the pitch angle between the k vector and B field (see Figj2fe/i). 
Accordingly, for CR with Larmor radius <C L, the injection scale of the 
turbulence, we obtain the pitch angle diffusion coefficient (Yan &; Lazarian 2007, 
henceforth YL07): 

dri^Jkmax{v)/L'n{'i- - V^) (2) 

where kmaxi't) is the cut-off wave number at ?], = is the cosine of wave pitch 
angle. Given D^^, the mean free path can be obtained by Imfp = 3/4 dfi "^^^^ ^ 



D 



v{l - /u"')^-^exp 



2^ 



3. Implications for Galactic CR transport and solar flare 

The scattering by fast modes is influenced by the medium properties as fast 
modes are usually subjected to damping, which varies from place to place. YL04 
has shown that the scattering properties for different ISM phases are different 
from each other. With the scattering at 90° known, we can make quantitative 
predictions for the mean free path of CRs in different phases. 
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Figure 3. Timescales involved in the acceleration in solar corona. The solid 
lines refer to the acceleration timescale. The dashed lines represent energy 
loss time scale. The dashdot lines give the escaping timescales. 



Galactic halo is in collisionless regime. The damping curve is given in 
Figl2|e/t. As we see, the damping decreases with the wave pitch angle 9. The 
mean free paths for moderate energy CRs are almost constant. This is due to 
the fact that gyroresonance changes marginally with the CR energy (see YL07 
for details). For higher energy CRs, the influence of damping saturates, and 
mean free path begins increasing with energy. In warm ionized medium (WIM) 
both viscous damping and collisionless damping are important (Figl2|e/t). We 
see that at low energies (< lOOGeV), the CR's mean free path decreases with the 
energy. This is because of the influence of viscous damping on gyroresonance. 

Since stochastic acceleration is intimately related to the scattering process, 
we apphed our result to acceleration for solar flares (Yan, Lazarian & Petrosian 
2007). We found that fast modes can provide efficient acceleration for both pro- 
tons and electrons (see Figj3]). In the case of protons, confinement of the particles 
can be achieved by fast modes as well, which is essential for the acceleration. 



4. Perpendicular transport 

In the three-dimensional turbulence, field lines are diverging away due to shear- 
ing by Alfven modes. Since the Larmor radii of CRs are much larger than the 
minimum scale of eddies /_L,mmj CRs will not be confined to the same field line. 
The cross-field transport thus results from the deviations of field lines as well as 
field line random walk on large scales (> min [L/M^, L]). 

High Ma turbulence: In high Ma turbulence, magnetic field becomes dy- 
namically important on the scale Ia = L/M^ (Lazarian 2006). If Ay <C Ia, the 
stiffness of B field is negligible so that D± = Dy. If Ay S> Ia, the diffusion is 
controlled by the straightness of the field lines, and D± = 1/31av- 

Low Ma turbulence: In this situation the steps in perpendicular direction 
are of LMl length. To diffuse over a distance R with random walk of LM^^ 
one requires {R/ LMj^)"^ steps. Therefore the perpendicular diffusion coefficient 
is /{R^ /D\\M\) = D\\M\. In the case Ay > L, the time of the individual 
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step is L/v\\, then D± = l/3LvM^. This is similar to the case discussed in the 

FLRW model (Jokipii 1966). However, we obtain the dependence of instead 
of their M\ scaling. 

Subdiffusion happens when particles are restricted to the magnetic field lines 
and the perpendicular transport is solely due to the random walk of field lines 
(see Kota &; Jokipii 2000). We do not have it in our earlier models. However, we 
may have sufficiently strong slab mode and it can cause subdiffusion. It is a very 
restrictive case. Consider a slab component. This component causes random 
walk of field lines in the plane perpendicular to the mean magnetic field direction. 
For a power-law spectrum E{k) oc {klsiab)~'^-i the critical scale up to which this 
component may induce random walk in the presence of Alfvenic field wandering 
is Zcrit = iklslab)~"^^ LlsiabM^^ Ag, where Ag = 5Bs/Bq is the dimensionless 
amplitude of the slab modes on scale Igiab- The dominant wavelength can be 
determined by maximizing 5bkk~^. Naturally, the constrain for this is that 
^crit- Otherwise no random walk has a meaning. Combining it with 

the expression of Zcriti we get a < acrit = 2. Detailed derivation provides 

1 

D± = AJslabD^ . 



5. Summary 

The cosmic ray transport is modified in accordance with the progress of un- 
derstanding of MHD turbulence. Implementing the tested model of turbulence, 
we found that fast modes dominate CR scattering if the turbulence energy is 
injected on large scales. On the other hand, the insights we gained through 
analytical studies give us guidance for designing future numerical experiments 
to test CR transport theories. 

Damping of fast modes leads to more dependencies of CR transport on the 
medium properties, different from the scattering by Alfven modes. Transient 
time damping (TTD) provides an important means of cosmic ray transport. Un- 
like gyroresonance, nonlinear broadening should be accounted. The 90° degree 
scattering can be realized by the broadened TTD, which ensures finite mean free 
paths for CRs. Gyroresonance with fast modes prevent streaming of moderate 
to high energy CRs, which cannot be confined by streaming instability. 

Perpendicular transport is diffusive, controlled by the diverging of field lines 
and field line random walk. The perpendicular diffusion depends on the CR's 
parallel mean free path and the strength of magnetic perturbation. Subdiffusion 
only applies below some critical scale if slab waves exist. 

In terms of numerical simulations of CR propagation, our findings of the 
dominance of fast modes, which, according to numerical studies are much more 
isotropic, compared to Alfven modes, may suggest that the present-day codes 
can capture important aspects of CR physics, provided that the damping of fast 
modes is properly taken into account. 
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